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ABSTRACT 

The ROSAT-ESO Flux Limited X-ray (REFLEX) galaxy cluster survey and the 2dF 
Galaxy Redshift Survey (2dFGRS) respectively comprise the largest, homogeneous X- 
ray selected cluster catalogue and completed galaxy redshift survey. In this work we 
combine these two outstanding datasets in order to study the effect of the large-scale 
cluster environment, as traced by X-ray luminosity, on the properties of the cluster 
member galaxies. We measure the Lx — cr^ relation from the correlated dataset and 
find it to be consistent with recent results found in the literature. Using a sample of 19 
clusters with Lx ^ 0.36 x 10^'* erg s~^ in the (0.1-2.4 keV) band, and 49 clusters with 
lower X-ray luminosity, we find that the fraction of early spectral type (77 ^ —1-4), 
passively-evolving galaxies is significantly higher in the high-Lx sample within i?2oo- 
We extend the investigation to include composite 6j cluster luminosity functions, and 
find that the characteristic magnitude of the Schechter-function fit to the early-type 
luminosity function is fainter for the high-Lx sample compared to the low-Lx sample 
(AM* — 0.58 ± 0.14). This seems to be driven by a deficit of such galaxies with 
Mjjj —21. In contrast, we find no significant differences between the luminosity 
functions of star-forming, late-type galaxies. We believe these results are consistent 
with a scenario in which the high-Lx clusters are more dynamically evolved systems 
than the low-Lx clusters. 

Key words: X-rays: galaxies: clusters - galaxies: evolution - galaxies: luminosity 
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1 INTRODUCTION 

It has been known since the earhest observations of rich clus- 
ters (see e.g., Abell 1965) that galaxies which inhabit these 
high-density regions are quite distinct from the general field 
galaxy population. The cluster population is dominated by 
galaxies with early-type morphologies, primarily ellipticals 
and SOs, which have colours and spectral types consistent 
with their undergoing passive evolution. In contrast, the field 
population is dominated by actively star-forming galaxies 
generally of late morphological types, such as spirals and 
irregulars. Explaining the segregation of galaxies into these 
two broad classes is one of the outstanding problems of ex- 
tragalactic astronomy and cosmology. 

Early work suggested that elliptical galaxies may have 
an intrinsically difi^erent formation process to disc galaxies 
(e.g., Sandage et al. 1970; Gott & Thuan 1976). However, 
since the general acceptance of the hierarchical process as 
the preferred model of structure formation (see e.g., Kauff- 
mann et al. 1993; Baugh et al. 1996; Cole et al. 2000), much 
attention has been focused on mechanisms that could trans- 
form star-forming, late-type galaxies into passively evolving, 
early-type systems in dense environments. Some examples of 
suggested evolutionary processes include: galaxy mergers - 
the results of numerical simulations suggest that the collision 
and merging of two equal mass disc galaxies can produce a 
product with properties typical of elliptical galaxies (e.g., 
Barnes 1992); ram-pressure stripping of gas from galaxies 
falling into clusters (Gunn & Gott 1972); and galajcy harass- 
ment, where tidal forces strip disc galaxies to make dwarf 
spheroidal systems (Moore et al. 1999). 

Modern large surveys such as the 2dF Galaxy Redshift 
Survey (2dFGRS, CoUess et al. 2001) and the Sloan Digital 
Sky Survey (SDSS, York et al. 2000), provide excellent qual- 
ity data and have allowed the dependence of galaxy proper- 
ties upon environment in the local universe to be examined 
in some depth in recent years. Lewis et al. (2002) studied the 
environmental dependence of the galaxy star formation rate 
near clusters in the 2dFGRS and found that it converged to 
the field value beyond ~ 3 times the virial radius, indicating 
that relatively small increments in local density at the out- 
skirts of clusters lead to a decrease in star formation rates. 
They found no dependence of the average star formation 
rate upon cluster velocity dispersion, a tracer of mass and 
hence a measure of the global environment. Gomez et al. 
(2003) used the Early Data Release of the SDSS to study 
the galaxy star formation rate as functions of local galaxy 
density and clustercentric radius. They found that the star 
formation rate is strongly correlated with the local galaxy 
density. Similarly to Lewis et al. (2002), Gomez et al. (2003) 
also found that the star formation rate converges to the field 
value at clustercentric distances of 3-4 virial radii. Balogh 
et al. (2004) extended these studies of 'galaxy ecology' to 
groups, using both 2dFGRS and SDSS data. They found 
little evidence that the distribution of the Ha line strength, 
an indicator of star formation rate, depends strongly upon 
environment amongst the actively star-forming galaxy pop- 
ulation. Similarly to both Lewis et al. (2002) and Gomez 
et al. (2003), they also found that the fraction of galaxies 
with significant ongoing star formation decreased steadily 
with increasing density, and interpreted the persistence of 
this correlation at low densities as indicating that ram pres- 



sure stripping is not the only mechanism responsible for the 
truncation of star formation in galaxies. In addition, Balogh 
et al. (2004) found little dependence of the star formation 
rate upon group velocity dispersion. 

The luminosity function (LF) of galaxies is a probabil- 
ity distribution over absolute magnitude, and by measuring 
the shape of the distribution in a number of environments, 
one is able to obtain clues as to which processes are impor- 
tant in shaping galaxy evolution. De Propris et al. (2003, 
hereafter DP03), constructed bj composite luminosity func- 
tions (LFs) of rich 2dFGRS clusters and found that the LFs 
of early spectral type galaxies have brighter characteristic 
magnitudes M* and steeper faint end slojics a than the 2dF- 
GRS field LFs of Madgwick et al. (2002). In contrast, the 
LFs of star-forming galaxies were found to be essentially 
identical. DP03 also constructed LFs for a variety of clus- 
ter subsamples, divided by several environmental variables 
including richness, velocity dispersion, substructure, and B- 
M (Bautz & Morgan 1970) type, but found no evidence of 
differences in the derived LF parameters. However, one mea- 
sure of the large scale cluster environment that DP03 were 
unable to test was cluster X-ray emission, which is produced 
from the thermal bremsstrahlung radiation of the hot intra- 
cluster gas (see e.g., Sarazin 1986), and is perhaps a better 
tracer of cluster mass than the velocity dispersion. There 
have of course been many other studies of the LF in differ- 
ent environments in recent years; we refer the reader to the 
references of DP03 for details of these. 

In this paper we seek to build upon the previous work 
described above and investigate the dependence of cluster 
galaxy populations upon the X-ray luminosity of the host 
cluster, by utilising the ROSAT-ESO Flux Limited X-ray 
galaxy cluster survey (REFLEX) catalogue (Boliringer et al. 
2004) in conjunction with galaxy data taken from the 2dF- 
GRS. This subject is also being explored by the RASS- 
SDSS galaxy cluster survey (Popcsso et al. 2004, 2005a, b,c), 
who have examined cluster scaling relations and galaxy lu- 
minosity functions for a sample of clusters selected from 
ROSAT X-ray observations, with galaxy photometry and 
spectroscopy drawn from the SDSS. 

The outline of this paper is as follows. In the next sec- 
tion, we briefly describe the two datasets used in this survey, 
2dFGRS and REFLEX. In Section 3, we describe our ini- 
tial method of correlating REFLEX and 2dFGRS data, and 
measure the relation between X-ray luminosity and cluster 
velocity dispersion. In Section 4, we supplement the RE- 
FLEX cluster sample with rich clusters taken from the cat- 
alogue compiled by De Propris et al. (2002, hereafter DP02), 
and divide the combined cluster sample in X-ray luminosity. 
We describe the construction of composite cluster luminosity 
functions in Section 5 for low- and liigli-Lx cluster subsam- 
ples, to complement the work of DP03. Finally, we end with 
a discussion of the results in Section 6 and our conclusions 
in Section 7. 

Throughout we assume a concordance cosmology of 
Qo = 0.3, f^A = 0.7, and Ho = 70 km s"^ Mpc"^ 

2 DATA 

Since the details of the 2dFGRS and REFLEX surveys have 
been covered in considerable depth elsewhere (see in par- 
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Figure 1. The bolometric Lx — Cr relation for the sample of 39 REFLEX clusters (dashed line). Only cluster members within of 
the cluster centroid were used in velocity dispersion estimation. 



ticular Colless et al. (2001); Norberg et al. (2002) for the 
2dFGRS, Bohringer et al. (2001); Collins et al. (2000) for 
REFLEX), here we only briefly review the gross character- 
istics of the datasets used in this paper. 



2.1 The 2dF Galaxy Redshift Survey 

The 2dFGRS is the largest completed spectroscopic survey 
of galaxies to date, covering approximately 1500 deg^ of the 
APM galaxy survey region (Dalton ct al. 1997), and contains 
spectroscopic information on approximately 250,000 galax- 
ies down to a magnitude limit of 6j — 19.45. The survey area 
consists primarily of two strips: the SGP strip centred close 
to the South Galactic Pole, covering an area of 80° x 15°, 
and the NGP strip in the Northern Galactic hemisphere that 
covers an area of 75° x 10° . In addition there are 99 widely 
scattered random 2° fields. The median redshift of the galax- 
ies in the survey is z = 0.11, and up to this redshift the 
survey is ~ 90 per cent complete at the stated 6j magnitude 
limit. 

In this paper we restrict ourselves to spectra with a 
2dFGRS quality fiag Q ^ 3. These galaxies have measured 

redshifts with an rms uncertainty of 85 km s~^ . In our exam- 
ination of galaxy populations, we make use of the 2dFGRS r) 
parameter defined by Madgwick et al. (2003). The value of "q 
is correlated with the present- to past-averaged star forma- 
tion rate, and so passively-evolving galaxies have low values 
of rj, and the opposite is true for galaxies that have formed 
a significant fraction of their stars more recently. In this pa- 
per we take a simple approach to classifying galaxies using 
this parameter, adopting r; ^ —1.4 for early- type, passively 
evolving galaxies (spectral type 1 of Madgwick et al. 2002), 
and J) > —1.4 for late- type, actively star- forming galaxies 
(encompassing spectral types 2-4 of Madgwick et al. 2002). 
Throughout this paper we adopt the following iiT-corrections 
for each population: 



K^^riy = 2.6z + 4.3«^ , and (1) 

KLte = l-3z + 2.0z\ (2) 

where -ff^ariy is the correction found for spectral type 1 
galaxies by Madgwick et al. (2002) and -ft'£ate is that found 
for spectral type 3, an appropriate average to use for the 
late-types. The rms error on the 2dFGRS 6j-band photom- 
etry is 0.15 mag. 

2.2 The ROSAT-ESO Flux Limited X-ray galaxy 
cluster survey 

REFLEX is the largest homogeneous catalogue of X-ray se- 
lected galaxy clusters assembled to date. REFLEX cluster 
candidates were selected from the region 6 < -|-2.5° of the 
second processing of the ROSAT All Sky Survey (RASS II, 
Voges et al. 1999) , avoiding the crowded stellar fields of the 
Magellanic clouds and the region ±20° either side of the 
Galactic plane. REFLEX therefore covers the entire 2dF- 
GRS region down to a nominal flux limit of 3 x 10~^^ erg 
s~^ cm~'^ in the (0.1-2.4 keV) energy band. The final cat- 
alogue (Bohringer et al. 2004) contains 447 clusters with 
measured redshifts, and is > 90 per cent complete. 

The mean fractional error on the measured REFLEX 
fluxes in the (0.1-2.4 keV) band is 16.7 per cent. Throughout 
this paper we quote REFLEX luminosities corrected for flux 
lost outside the detection aperture as described by Bohringer 
et al. (2000). These can be found in Table 6 of Bohringer 
et al. (2004). 

3 CORRELATING REFLEX AND THE 2dFGRS 

3.1 Membership selection within i?v 

In our initial attempt at cluster membership determination, 
we selected 2dFGRS galaxies located within the virial radius 
i?v of each REFLEX cluster. We calculated Rv by adopting 
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Table 1. Recent measurements of the L^j^^^^ — ar relation, in the form L-j^^^^ = 10™ x (jT.(km s ^) erg s ^. ATeiugters is the number of 
clusters in each sample. 



Reference 



Comments 



White, Jones & Forman (1997) 


5.36±0.16 


og 3-H0.13 
9 


14 


Cooling flow clusters removed. 


Mahdavi & Geller (2001) 




31.81°:^ 


280 




Girardi & Mezzetti (2001) 


4.4ll;g 


29.4t3;0 


51 


Multicomponent clusters removed. 


Xue & Wu (2000) 


5.30±0.21 


28.32±0.61 


197 




Ortiz-Gil ct al. (2004) (a) 


4.1±0.3 


32.72±0.08 


171 




Ortiz-Gil et al. (2004) (b) 


4.2±0.4 


32.41±0.10 


123 


Multicomponent clusters removed. 


This work 


4.8±0.7 


30.6±2.1 


39 





a self-similar model, taken from Arnaud et al. (2002) and 
references therein. In the simplest self-similar assumption, 
Rv depends only on the cluster X-ray temperature and a 
fixed density contrast with respect to the critical density of 
the universe at redshift z. We converted the (0.1-2.4 keV) 
REFLEX X-ray luminosities to bolometric luminosities us- 
ing Table 5 in the REFLEX catalogue paper (Bohringcr 
et al. 2004), the X-ray temperature of each cluster being 
estimated using the I'x(o.i-2.4keV) — T relation of Maxke- 
vitch (1998), uncorrected for cooling flows since this was 
the relation used to calculate the REFLEX (0.1 2.4 keV) 
band X-ray luminosities. The bolometric Lx — T relation of 
Markevitch (1998) was then applied to estimate the X-ray 
temperature used in the calculation of (however, in prac- 
tice this extra step has a negligible effect on the estimated 
value of R^). 

Membership along the line of sight was determined in 
an iterative fashion, using the biweight scale estimator for 
velocity dispersion cr,. recommended by Beers et al. (1990). 
Initially, ar was calculated for galaxies within ±2000 km s~^ 
of each REFLEX cluster redshift. On subsequent iterations, 
a conservative 3ar clipping was applied, and galaxies with 
line of sight velocities outside of this range with respect to 
the cluster redshift were discarded. The procedure was found 
to converge within a few iterations, and at the end of the 
process clusters with less than 15 members were discarded, 
since the calculated velocity dispersions in these cases were 
unlikely to reflect the true values. The uncertainty in (Jr 
was estimated via bootstrap resampling, and was found to 
be typically ~ 80 km s^^. In total, we extracted 39 RE- 
FLEX clusters with 15 members from the 2dFGRS using 
this method. Typically, each REFLEX cluster contains ~ 80 
2dFGRS member galaxies. 



3.2 The ix - cTr relation from 2dF-REFLEX data 

Since the X-ray emitting gas and galaxies that make up a 
cluster share a common potential well, a relationship exists 
between Lx and ffr- Self-similar models, where it is assumed 
that the evolution of clusters is solely due to their collapse 
under gravity, predict that I/Xboi °^ '^r (see e.g. Quintana 
& Melnick 1982). Some authors have obtained significantly 
steeper values (e.g.. White, Jones & Forman 1997; Xue & Wu 
2000), and these results have been interpreted as evidence 
that feedback (i.e. non-gravitational heating) processes are 
important in clusters. 

The most relevant recent work in our case is that of 
Ortiz-Gil et al. (2004), who measured the Lx — <Tr relation 



for the REFLEX team using a combination of literature and 
optical follow-up data. We measured the relation using the 
39 cluster 2dFGRS-REFLEX sample, to investigate whether 
the homogeneity of the combined dataset could overcome 
the small sample size and provide competitive constraints 
on this relation. The results for bolometric luminosities are 
shown in Fig. 1. We used the BCES (bivariate correlated er- 
rors and intrinsic scatter) bisector method of Akritas & Ber- 
shady (1996) to obtain the coefficient and power-law slope 
estimates of the relation. This fitting technique takes into 
account errors in both variables and intrinsic scatter. The 
following results were obtained: 



in30.6±2.1 /, _i^4.8±0.7 _i 

= 10 ar (km s ) erg s 



and 



(3) 



Lx(o.i-2.4 kov) = 10 (Tr(km s ) erg s . (4) 

The uncertainties quoted are the la errors determined 
by bootstrap resampling, and their large size places only 
weak constraints on the relation these results arc con- 
sistent with measurements that favour feedback, but also 
do not rule out simple, self-similar evolution (see Table 1). 
This is true for many of the results listed in Table 1 with 
the exception of White et al. (1997) and Xue & Wu (2000), 
which differ significantly from the self-similar prediction for 
reasons that are not clear. Table 1 suggests that the large 
amount of intrinsic scatter in this relation can only be over- 
come by using very large sample sizes and/or cleaning the 
cluster sample of objects containing cooling flows. 



4 SUPPLEMENTING THE REFLEX CLUSTER 
SAMPLE 

4.1 Membership selection within R200 

Of the 39 REFLEX clusters extracted from the 2dFGRS 
using the method described in Section 3, only 23 of these 
clusters are located within the redshift range z < 0.11, for 
which the 2dFGRS is nearly complete, and have an over- 
all 2dFGRS completeness (i.e. the fraction of galaxies from 
the parent APM galaxy catalogue with measured rcdshifts) 
at the cluster centroid of > 70 per cent. We found that this 
sample was too small to practically subdivide in X-ray lumi- 
nosity in order to study the effect of Lx on the cluster galaxy 
luminosity function. De Propris et al. (2002, hereafter DP02) 
have previously studied Abell, EDCC and APMCC clusters 
within the 2dFGRS region, and within z < 0.11 their cata^ 
logue should represent am essentially complete sample of rich 
clusters. REFLEX nominally contains all the clusters in the 
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Table 2. The high-Lx cluster sample. Nmem is the number of 2dFGRS cluster member galaxies with redshift quality flag Q ^ 3 used 
in the estimation of <Tr- 



High-Lx 
0.36 X lO'*'* erg s'^) 


RA 
(Deg. 2000) 


Dec. 
(Deg. 2000) 


z 


(km s-1) 


^200 

(Mpc) 


-'^X(0.1-2.4 keV) 

(xlO** erg s-i) 


^ mem 


Completeness 


A0954 


153.4370 


-0.1085 


0.0950 


780±80 


1.9 


0.764±0.114 


48 


0.86 


A0957 


153 4180 


-0.9144 


0.0451 


730-1-50 


1.8 


n 434-1-0 056 


85 


0.85 


A1650 


194.6710 


-1.7569 


0.0841 


700±60 


1.7 


3.873±0.236 


97 


0.84 


A1651 


194 8390 


-4.1948 


0.0842 


940±70 


2.2 


4 289-1-0 257 


115 


0.80 


A1663 


195.7110 


-2.5062 


0.0829 


810±70 


1.9 


0.815±0.178 


88 


0.84 


A1750 


202.7080 


-1.8728 


0.0858 


920±50 


2.2 


2.291±0.266 


116 


0.82 


A2734 


2.8364 


-28.8551 


0.0612 


810±60 


2.0 


1.197±0.108 


108 


0.82 


A2811 


10.5362 


-28.5358 


0.1078 


950±70 


2.2 


3.030±0.297 


94 


0.93 


A3027 


37.6814 


-33.0987 


0.0774 


910±70 


2.2 


0.449±0.070 


90 


0.75 


A3094 


47.8540 


-26.8998 


0.0683 


710±60 


1.7 


0.362±0.182 


90 


0.90 


A3880 


336.9690 


-30.5699 


0.0577 


820±60 


2.0 


0.939±0.077 


105 


0.96 


A4038 


356.9300 


-28.1414 


0.0303 


910±40 


2.2 


1.127±0.043 


155 


0.91 


S1136 


354.0710 


-31.6103 


0.0620 


650±70 


1.6 


0.549±0.113 


43 


0.89 


S0084 


12.3502 


-29.5244 


0.1087 


820±70 


1.9 


1.566±0.251 


52 


0.93 


S0041 


6.3849 


-33.0472 


0.0494 


580±50 


1.4 


0.537±0.052 


70 


0.81 


RXCJ1326.2+13 


201.5740 


0.2257 


0.0827 


560±110 


1.3 


1.001±0.117 


31 


0.81 


RXCJ1309.2-136 


197.3210 


-1.6126 


0.0831 


560±70 


1.3 


1.030±0.156 


30 


0.88 


RXCJ0229.3-3332 


37.3430 


-33.5378 


0.0774 


760±60 


1.8 


0.606±0.078 


58 


0.73 


RXCJ0225. 1-2928 


36.2939 


-29.4740 


0.0607 


550±60 


1.3 


0.434±0.102 


53 


0.73 



surveyed region above a flux limit of 3 x 10~^^ erg cm~^ 
in the (0.1-2.4 keV) energy band. We therefore decided to 
combine the non- REFLEX clusters catalogued in DP02 with 
REFLEX, and use the REFLEX survey flux limit to divide 
the combined cluster sample in terms of X-ray luminosity. 

To determine the cluster membership of the whole sam- 
ple in a consistent way, we chose to replace the virial radius 
as used in Section 3 with R200, defined as the radius within 
which the mean interior density is 200 times the critical 
density. This is straightforward to calculate by the equation 
(from Carlberg et al. 1997) 



R200 — 



mH{z)' 



(5) 



where H{z) is the Hubble constant at redshift z. 

Since R20Q is defined by the measured cluster velocity 
dispersion, it was necessary to calculate both these quanti- 
ties iteratively, in step with one another. Therefore, on the 
first iteration the velocity dispersion was calculated using 
all galaxies within a 1-Mpc radius of the catalogued cluster 
centre. The initial velocity cut applied around the cluster 
redshift was ±2000 km s~^ as before. After each iteration, 
R200 was calculated using equation (5) and the 3ar velocity 
clipping was applied. This process again converged rapidly, 
within 5 iterations. 

As before, clusters with less than 15 members at the end 
of the process were discarded. This left 36 REFLEX clus- 
ters, with measured velocity dispersions in agreement within 
the uncertainties of the values obtained using the selection 
method described in Section 3 (i.e. 2 clusters with few mem- 
bers were lost compared to using _Rv). R200 is often taken to 
be equivalent to the virial radius, and it was found from the 
REFLEX sample that typically ii200 ~ (0.9±0.2)iiv, where 
the quoted uncertainty is the standard deviation. Galax- 
ies within 7i2oo are therefore likely to be gravitationally 
bound to the cluster. 146 out of the 284 clusters catalogued 
by DP02 were recovered by the method. A comparison of 



the member numbers of the Abell clusters obtained by our 
method with Table 1 of DP02, who used a different method 
for determining cluster membership based on a 'gappcr' al- 
gorithm (described in Beers et al. 1990), revealed that the 
clusters lost in the process were located typically at ^ > 0.11, 
and had less than the specified minimum 15 cluster mem- 
bers. We found generally good agreement with the velocity 
dispersions calculated using the biweight scale estimator in 
comparison to the DP02 values - 70 per cent of the 146 
clusters recovered had velocity dispersions within < Icr of 
agreement. 



4.2 Defining the cluster samples 

In the following, we consider only clusters with an overall 
2dFGRS redshift completeness > 70 per cent and z < 0.11 
- since within the survey median redshift the 2dFGRS has 
a very high redshift completeness up to the survey magni- 
tude limit. These criteria result in a sample of 23 REFLEX 
clusters and 94 DP02 non-REFLEX clusters. 

To divide the whole cluster sample in terms of X-ray 
luminosity, we make use of the fact that REFLEX is a flux 
limited survey, assuming that the REFLEX catalogue con- 
tains all clusters above the survey flux-limit. We estimated 
a maximum X-ray luminosity for each non-REFLEX cluster 
using the nominal flux-limit of the REFLEX cluster cata- 
logue (3 X 10"^^ erg s"^ cm~^ in the (0.1-2.4 keV) energy 
band) . We chose to divide the combined cluster sample into 
two subsamples, defined by 



ix(o.i-2.4 kev) < 0.36 X 10** erg s ^ and 

Lx{i).i-2A keV) > 0.36 X lO'^'crg s"^ , 



(6) 
(7) 



for the low- and high-Lx samples respectively. The Lx value 
around which to divide the sample was chosen to maximise 
the size of the low-Lx cluster sample, while ensuring ade- 
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Table 3. The low-Lx cluster sample. Nmem is the number of 2dFGRS cluster member galaxies with redshift quality flag Q ^ 3 used in 
the estimation of ffr- Clusters marked with asterisks (*) were identified by De Propris et al. (2002) as being located along the same line 

of siglit at (lifForoiit rodshifts. 
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quatc coverage of the galaxy magnitude range for the con- 
struction of composite LFs for the high-Lx sample. 

The high-Lx sample consists of 19 REFLEX clusters, 
while the low-Lx sample is composed of 49 clusters, 4 from 
REFLEX with the remainder coming from the catalogue of 
DP02. Because of the way in which wc have defined these 
samples, the low-Lx clusters can be thought of as essentially 
an optically-selected sample that has been cleaned of clus- 
ters that could potentially have X-ray luminosities above 



the division we have made. In the following analysis, we 
are therefore investigating the effect that selection of clus- 
ters based upon the observation of hot. X-ray emitting gas 
has on the member galaxy populations. We assume that the 
DP02 clusters are bound systems, since they were selected 
from a redshift survey and should not be artifacts of projec- 
tion effects. The properties of the high- and low-Lx clusters 
are summarised in Tables 2 and 3 respectively. 
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Figure 2. (a) Line of sight velocity dispersion distributions for the low- and high-Lx cluster samples, (b) Redshift distributions. 
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Figure 3. Mean fraction of carly-typc (r; —1-4) galaxies versus fraction of R200 for the low- and high-Lx cluster samples. A volume- 
limited galaxy sample, defined by M;,^ < —19, z < 0.11, was used. The error bars are the standard error on the mean in each bin. 



4.3 Properties of the cluster samples 

Having divided the sample in terms of X-ray luminosity on 
the basis of several assumptions, it is useful to check the 
velocity dispersion distributions of the resulting low- and 
high-Lx samples, because Lx and ar are correlated. Panel 
(a) of Fig. 2 shows that the majority of the low-Lx sample 
have CTr < 600 km s~^, whereas in contrast most of the 
clusters in the high-Lx sample have ar > 700 km s~^. The 
median values of the low- and high-Lx samples were 510 
km s^^ and 780 km s^^ respectively. This is consistent with 
the high-Lx sample on average containing the higher mass 
systems, as expected. Panel (b) of Fig. 2 shows the redshift 
distributions of both samples, and indicates that the average 
redshift of the high-Lx clusters is slightly higher than that 
for the low-Lx clusters. This is primarily due to the way in 
which we have defined our samples, but it is also in part due 
to the fact that the high-Lx clusters are intrinsically rarer 
systems. The median redshifts were 0.056 and 0.077 for the 
low- and high-Lx samples respectively. 

Under the assumption that the redshift completeness of 
the 2dFGRS is not a function of galaxy spectral type within 
z < 0.11, we are immediately able to examine the mix of 
spectral types in each cluster sample. In Fig. 3, we plot the 
mean fraction of passively evolving, early- type (77 ^ —1.4) 
galaxies in composite low- and high-Lx clusters in 20 per 



cent radial bins of L?200- Since our samples have different 
median redshifts, we remove the (Icpciidciicc upon magni- 
tude by using a volume-limited galaxy sample defined by 
Mbj < -19, z < 0.11. The error bars in Fig. 3 are the 
standard errors on the mean galaxy fraction in each bin. 
Clearly, the fraction of passively evolving, early-type galax- 
ies is higher at all radii out to 7?20() in the high-Lx sample 
compared to the low-Lx sample. In both cases, the frac- 
tion of early-type galaxies falls off smoothly with increasing 
radial distance from the cluster centre. The overall mean 
early-type galaxy fraction within i?200 is 0.76±0.02 for the 
high-Lx sample compared to 0.64±0.02 for the low-Lx sam- 
ple. 



5 COMPOSITE LUMINOSITY FUNCTIONS 
5.1 Construction 

We decided to examine the implications of the different mix 
of spectral types found in the low- and high-Lx clusters on 
the composite cluster luminosity functions by following the 
method of CoUess (1989). The summation carried out for 
each absolute magnitude bin of the composite LF was 

A^co 



£Vco \ " Nij_ 
m, 2^ Nio ' 



(8) 
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Figure 4. Composite luminosity functions for all spectral types for: (a) the high-Lx sample and (b) the low-Lx sample. The solid line 
in each plot is the best fitting Scliechtcr function. For comparison, we plot the renormalised fit to the liigh-Lx sample LF as the dashed 
line in panel (b). Quoted independent errors on each fitted parameter are given at the la level, and were determined using a bootstrap 
resaniijling technique, (c) The deviation of each luminosity function from the best-fitting Schechter function fit to the low-Lx sample 
data, in units of the error on each point, (d) 1, 2, 3a error contours for joint luminosity function parameter estimates of the low- and 
high-Lx samples. 



Here Ncj is the number of galaxies in the jth bin of the 
composite LF, Nij is the number of gala^xies contributing 
to the jth. bin of the ith. cluster LF, and rrij is the number 
of clusters contributing to the jth bin of the composite LF. 
A^io is a normalisation parameter (effectively a richness) for 
the ith cluster LF, in our case calculated from the complete- 
ness corrected number of galaxies with Mtj < —19, which 
defines a volume-limited sample of 2dFGRS galaxies within 
z < 0.11. N'co is simply the sum of all the individual cluster 
normalisations: 



sured by the 2dFGRS and ii'-corrected them accordingly 
using equations (1) and (2). 

The errors in were calculated according to: 



SNcj = 



^ V Mo 



1/2 



(10) 



(9) 



In contrast to DP03, who used the parent 2dFGRS pho- 
tometric catalogue (i.e. the APM galaxy survey) directly to 
correct for incompleteness, we corrected Nij for overall and 
magnitude-dependent incompleteness using the 2dF survey 
masks and software as described in Appendix A of Norberg 
et al. (2002). We found that the average completeness correc- 
tions across all bins were 15 per cent for the high-Lx sample 
and 11 per cent for the low Lx sample. In constructing all of 
our LFs, we used the extinction corrected magnitudes mear 



5.2 Results 

We constructed LFs for galaxies of all, early, and late spec- 
tral types separately, for both the low- and high-Lx cluster 
samples. Each LF was divided into 10 bins over the magni- 
tude range —17.3 > Mi,j > —22.3. This ensured that at least 
three clusters contributed to each LF bin. The faint magni- 
tude limit that we adopt for our LFs is 0.5-mag brighter 
than that adopted by DP03. 

We fitted Schechter (1976) functions. 



nc{M)dM = kN*e 



k{M*-M) 



dM, 



(11) 



to each LF, where M is the absolute 6j magnitude, k = 
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Figure 5. Composite luminosity functions for early (r/ ^ —1.4) spectral types for: (a) the higli-Lx sample and (b) the low-Lx sample. 
The solid line in each plot is the best fitting Schcehtcr function. For comparison, we plot the renormalised fit to the high-Lx sample 
LF as the dashed line in panel (b). Quoted independent errors on each fitted parameter are given at the la level, and were determined 
using a bootstrap resampling technique, (c) The deviation of each luminosity function from the best-fitting Schechter function fit to the 
low-Lx sample data, in units of the error on each point, (d) 1, 2, 3<t error contours for joint luminosity function parameter estimates 
of the low- and high-Lx samples. 



ln(10)/2.5, M* is the characteristic magnitude of the turn- 
off at the bright-end and a is the faint-end slope. We per- 
formed two-parameter fits to each LF for M* and a using a 
minimisation technique. A'^* was fixed such that the value 
of the integral of the fitted function was equal to the total 
number of galaxies in the composite LF within the same 
magnitude range. We also tried fitting for the normalisation 
with N* held as a free parameter, but found that this made 
no significant difference to our results. 

The la errors on the individual Schechter function pa- 
rameters were estimated using a bootstrap resampling tech- 
nique, with the other parameter fixed to the value obtained 
in the minimisation. In each bootstrap sample, we forced 
the number of galaxies drawn from each cluster to be equal 
to the number of galaxies found in that cluster in the real 
sample. This ensured that the weighting of each cluster in 
the production of the composite LF was approximately the 
same for each bootstrap sample. We believe that the er- 
ror estimates obtained using this technique are more robust 
than X confidence intervals for individual parameters, be- 
cause by design the bootstrap method takes into account 
the effect of sample-to-sample variation. We found that the 
errors estimated using the bootstrap resampling technique 



were typically ~ 50 per cent larger than the equivalent Icr 
individual parameter confidence intervals. 

We now examine the results for each spectral binning 
in turn. 



5.2.1 All spectral types 

Fig. 4 shows the resulting composite LFs for galaxies of all 
spectral types in the low-Lx and high-Lx samples. Panel 

(a) of this figure shows that the Schechter function is a poor 
fit to the high-Lx sample LF at the bright-end, where there 
is a significant excess of galaxies above that predicted by 
the fitted model. The two brightest bins of this LF contain 
the Brightest Cluster Galaxies (BCGs) of several clusters, 
and we found that explicitly removing the BCGs from the 
sample made no significant difference (< Icr) to the derived 
LF parameters. Note that to make a fair comparison with 
the low-Lx LF we have fitted the high-Lx LF over the same 
magnitude range. However, fitting the high-Lx LF over the 
full range in absolute magnitude ( — 17.3 > Mb, > —23.3) 
does not significantly alter the values of the Schechter func- 
tion parameters. We obtain M* = -20.07 ± 0.06, a = 
—0.84 ± 0.05 when the very brightest galaxies are included. 
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Figure 6. Composite luminosity functions for late (rj > —1.4) spectral types for: (a) the liigli-Lx sample and (b) the low-Lx sample. 
The solid line in each plot is the best fitting Schechter function. For comparison, we plot the renormalised fit to the high-Lx sample 
LF as the dashed line in panel (b). Quoted independent errors on each fitted parameter are given at the Icr level, and were determined 
using a bootstrap resampling technique, (c) The deviation of each luminosity function from the best-fitting Schechter function fit to the 
low-Lx sample data, in units of the error on each point, (d) 1, 2, 3a error contours for joint luminosity function parameter estimates 
of the low- and high-Lx samples. 



In contrast, the Schechter function provides a good fit to the 
LF of the low-Lx clusters along the whole range in absolute 
magnitude, as can be seen from panel (b). 



In panel (c), we show the difference between the LF 
value at each point and the best-fitting Schechter frinction 
fit for the low-Lx sample (after renormalisation to the high- 
Lx LF). The vertical axis in this plot is in units of the error 
on each LF data point. We can see from this plot that there 
is a significant deficit of galaxies at Mi,j ~ — 21 in the high- 
Lx sample LF in comparison to the best-fitting Schechter 
function for the low-Lx sample. Panel (d) shows the joint 
1, 2, 3a error contours on the two-parameter fit. Some 
overlap of the confidence contours at the 2a level is visible 
in this plot. When we consider the uncertainties on the in- 
dividual parameter estimates, we find significant differences 
between the low- and high-Lx LFs. M* is fainter in the 
high-Lx sample than the low-Lx sample at 4.7(7 significance 
(AM* = 0.51±0.11). The faint-end slope is also shallower in 
the high-Lx sample at 3.3a significance (Aa = 0.23 ±0.04). 



5.2.2 Early spectral types 

Because both cluster samples are dominated by passively 
evolving galaxies, the results for the early-type LFs shown 
in Fig. 5 follow a similar pattern to that for the LFs of 
all spectral types. We see again that the Schechter function 
is a poor fit to the high-Lx sample LF at the bright-end, 
and that the low-Lx sample is well fitted by the Schechter 
function over its entire range. Again, a fit to the high-Lx LF 
over the full range in absolute magnitude (—17.3 > Mbj > 
—23.3) does not give significantly different results from those 
obtained for the range (—17.3 > Mbj > —22.3) - we obtain 
M* = -20.14 ± 0.08, a = -0.80 ± 0^06. M* and a are both 
found to be lower in the high-Lx clusters compared to the 
low-Lx sample. 

Inspection of panel (d) shows a small amount of overlap 
between the 2a confidence contours for the joint parameter 
estimates. If we consider the uncertainties on the individual 
parameters, the difference in the LF characteristic magni- 
tudes is significant at the Aa level (AM* = 0.58 ± 0.14). In 
contrast, the difference between the a values is only sig- 
nificant at the 2.3<T level (Aa = 0.18 ± 0.08). Panel (c) 
shows the main reason for the lower value of M* in the 
Schechter function fit to the high-Lx sample LF - there is a 
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Table 4. Comparison of derived LF parameters with the values of De Propris et al. (2003) and Madgwick et al. (2002), based on 2dFGRS 
data. Quoted errors on the individual parameter estimates for this work are Icr errors determined using a bootstrap resampling technique. 



Reference 


Spectral type 


A/T* 
IVl 


Ot 


Sample 


This work — low-Lx 


All 


-20.57±0.09 


-1.07±0.05 


49 clusters 


This work - liigli-Lx 


All 


-20.06±0.06 


-0.84±0.05 


19 clusters 


De Propris et al. (2003) 


All 


-20.84±0.07 


-1.28±0.03 


60 clusters 


Madgwick et al. (2002) 


All 


-20.56±0.04 


-1.19±0.01 


Field 


This work — low-Lx 


Early 


-20.70±0.12 


-0.97±0.06 


49 clusters 


This work - higli-Lx 


Early 


-20.12±0.08 


-0.79±0.06 


19 clusters 


Dc Propris ct al. (2003) 


Type 1 


-20.81±0.09 


-1.05±0.04 


60 clusters 


Madgwick et al. (2002) 


Type 1 


-20.35±0.05 


-0.54±0.02 


Field 


This work - low-Lx 


Late 


-20.07±0.14 


-1.09±0.08 


49 clusters 


This work - high-Lx 


Late 


-19.79±0.11 


-0.93±0.10 


19 clusters 


De Propris et al. (2003) 


Type 2 


-20.25±0.13 


-1.23±0.07 


60 clusters 


De Propris et al. (2003) 


Types 3-1-4 


-19.91±0.19 


-1.30±0.10 


60 clusters 


Madgwick et al. (2002) 


Type 2 


-20.30±0.03 


-0.99±0.01 


Field 


Madgwick et al. (2002) 


Type 3 


-19.94±0.04 


-1.24±0.02 


Field 


Madgwick et al. (2002) 


Type 4 


-19.92±0.05 


-1.50±0.03 


Field 


Madgwick et al. (2002) 


Late (weighted average) 


-20.11±0.04 


-1.16±0.02 


Field 



significant (~ 4<t) deficit of early spectral type galaxies with 
Miyj ~ —21 in comparison to the best Schechter function fit 
to the low-Z/x sample LF. 

5.2.3 Late spectral types 

In panels (a) and (c) of Fig. 6, we see that the high-Lx 
LF contains no galaxies with Af[,, < —21.3. Despite this, 
panel (d) shows that there are no significant differences be- 
tween the LF parameter values obtained for each sample 
when the confidence contours for the joint parameter fit are 
considered. The uncertainties in the individual LF parame- 
ters indicate likewise: AM* = 0.28 ±0.18, Aq = 0.16 ±0.13. 
Therefore the differences in the LF parameters found for 
galaxies of all spectral types axe caused by the different LFs 
of early- type, passively evolving galaxies. 



6 DISCUSSION 

6.1 Early-type galaxy fractions 

In Section 4.3 we established that in both the low- and high- 
Lx cluster samples the early spectral type galaxy fraction 
falls steadily with increasing distance from the cluster cen- 
tre (Fig. 3). Although the 2dFGRS spectral type parame- 
ter 77 is correlated with morphology with a large amount of 
scatter, this behaviour is reminiscent of the long-established 
galaxy morphology-local density relation in clusters (e.g. 
Dressier 1980). In Fig. 3, we also see that the fraction of 
passively evolving, early spectral type galaxies is higher in 
the high-Lx sample than in the low-Lx sample. Balogh et al. 
(2002) reported a result similar in character to that found in 
this study from a morphological study of low- and high-Lx 
clusters conducted with the Hubble Space Telescope. They 
found that within a fixed physical region, the mean fraction 
of disc-dominated galaxies is strongly dependent upon Lx. 
Although the epoch of this study was different {z ~ 0.25 
compared to 2: < 0.1 for our study), as are the observations 



and quantities measured, the fact that spectral type is corre- 
lated with morphology, albeit with large scatter (Madgwick 
et al. 2003), means that a comparison of our result with 
one of the conclusions of Balogh et al. (2002) is not un- 
reasonable. Balogh et al. (2002) interpreted their results as 
indicating that galaxy mergers have played a bigger role in 
the evolution of more massive clusters. 

Since both Lx and Oy should trace the cluster mass, and 
7) is correlated with the equivalent width of the Ha line, we 
can also mate a comparison with the results of Balogh et al. 
(2004), who examined the correlation between the fraction 
of galaxies with Hq equivalent width > 4A and local den- 
sity in the 2dFGRS and SDSS. Balogh et al. (2004) divided 
a sample of ~ 25, 000 group and cluster galaxies into two 
subsamples of low-cTr (200 < Ur < 400 km s~^) and high-Ur 
(500 < CTr < 1000 km s~^). They found marginal evidence 
that the fraction of star-forming galaxies was lower in the 
high-CTr clusters relative to the low-cr,- groups at fixed local 
density. This is similar to the result we see for the low- and 
high-Lx cluster samples studied in this paper, although we 
make no attempt to distinguish between the effects of varia- 
tion in local galaxy density and the large scale cluster envi- 
ronment (as traced by Lx). Lewis et al. (2002) also studied 
star-formation rates in 2dFGRS cluster galaxies as traced by 
the equivalent width of Ha. Lewis et al. (2002) considered 
17 clusters drawn from DP02 and divided their sample into 
clusters with Gr > 800 km s"'^ and ar < 800 km s~^ , with 10 
and 7 clusters in each respective subsample. They found no 
dependence of galaxy star-formation rates upon the velocity 
dispersion of the host cluster, in contrast to the results that 
we obtain for our low- and high-Lx cluster samples. This is 
perhaps due to the comparatively small size of the cluster 
sample used by Lewis et al. (2002). 



6.2 Luminosity functions 

We present a comparison of our derived LF parameters with 
the values obtained from 2dFGRS data by DP03 for clus- 
ters and Madgwick et al. (2002) for the field in Table 4, 
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Table 5. Luminosity function parameters derived for studies using SDSS g-band photometric data, transformed to Hq = 70 km s 
Mpc-i and the 2dFGRS bj band. 



Reference 


Spectral type 


M* 


a 


Sample 


Popesso et al. (2005a) 

Popcsso et al. (2005c) 
Blanton ct al. (2003) 


All 
All 
All 


-21.11±0.20 
-20.70±0.21 
-19.91±0.02 


-1.33±0.04 
-1.07±0.12 
-0.89±0.03 


69 clusters (members selected within 1.5 Mpc radius) 
69 clusters (members selected within ii20o) 

Field 


Popesso et al. (2005c) 


Early 


-20.31±0.16 


-0.69±0.10 


15 clusters (members selected within -R200) 


Popesso et al. (2005c) 


Late 


-21.65±0.40 


-1.80±0.04 


15 clusters (members selected within R200) 



divided by spectral type and transformed to Ho = 70 km 
s~^ Mpc~^. Here, spectral type 1 corresponds to our defini- 
tion of early-type galaxies and spectral types 2-4 collectively 
form our sample of late-type galaxies. In this table we also 
quote average values of M* and a for the late-type LFs of 
Madgwick et al. (2002), weighted by the respective normal- 
isation constants, to ease comparison with our late-type LF 
parameter values. 

The most surprising result of this work is the disagree- 
ment of the high-Lx sample early-type LF parameter val- 
ues with both the low-Lx sample and the results of DP03. 
Al* and a are both lower in the high-Lx, carly-typc LF 
in comparison to the corresponding DP03 early-typo LF 
at > 3(T significance. The characteristic magnitude of the 
high-I/x sample LF is also significantly lower than found in 
the field by Madgwick et al. (2002), although the faint-end 
slope is steeper. In contrast, both Schcchtcr function param- 
eter values for the early-type LF of the low-Lx sample are 
within ~ la of agreement with the results of DP03, with 
M* brighter than in the field, and steeper a. 

The late-type LF parameter values for the high-Lx sam- 
ple agree within < 3a of the typo 2 and types 3+4 LFs of 
DP03, and the weighted average of the field results of Madg- 
wick et al. (2002). There is therefore marginal disagreement 
- although note that there is no significant disagreement 
between the low- and high-Lx late-type LFs. The low-Lx 
sample late-typo LF parameters are well within < 2a of 
agreement with the DP03 results and the Madgwick et al. 
(2002) field results. These results are therefore consistent 
with a picture in which the late-type LF varies little with 
environment. 

We are also able to compare our results to similar stud- 
ies undertaken using SDSS g-band data, following a suitable 
transformation to the bj band as used by the 2dFGRS, i.e., 

bj = g + 0.15 + 0.13{g-r), (12) 

as used by De Propris et al. (2004). We present results from 
Popesso et al. (2005a, c) for RASS-SDSS galaxy clusters, and 
Blanton et al. (2003) for the field in Table 5. Note that we 
assume a (g—r) colour of 0.8 when converting the M* values 
using equation (12), a value we estimated from inspection 
of the colour distribution of early spectral type 2dFGRS 
galaxies, after transforming the 2dFGRS bj and Rf bands 
to Sloan g and r using the prescriptions given in De Propris 
et al. (2004). 

Popesso ct al. (2005a) computed LFs for a sample of 69 
X-ray selected clusters drawn from the RASS-SDSS survey. 
We quote the LF parameter values for the composite LF 
produced from galaxies within 1.5 Mpc clustercentric dis- 



tance in Table 5. We can see from comparison of the results 
presented in Section 5.2 that there is a significant disagree- 
ment between Popesso et al. (2005a) and the results that we 
obtain for both LF parameters for the high-Lx sample LF 
of all spectral types - in fact, we obtain a value of M* that 
is ~ 1 magnitude fainter. In contrast, there is mild disagree- 
ment (at the ~ 2.5a level) between the characteristic magni- 
tude of this LF and the value that we obtain for the low-Lx 
sample LF for galaxies of all spectral types. However, our 
low-Lx sample LF has a significantly fiatter faint-end slope 
than this LF. 

Popesso ct al. (2005c) constructed LFs for the same 
69 RASS-SDSS clusters using galaxies selected within J?2oo, 
rather than using a fixed metric aperture, and fitted a com- 
posite of two Schechter functions to each LF. Using this 
method, they obtained results consistent with those previ- 
ously obtained in Popesso et al. (2005a). It is appropriate 
to compare the LF values that we obtained in Section 5.2 
with the bright component of the Popesso et al. (2005c) 
double-Schechter function fits. We find in this case a much 
better agreement of our low-Lx sample LF for galaxies of 
all spectral types, with the values of M* and a in agree- 
ment within ~ la of each other. The agreement with the 
LF parameters we derived for the equivalent high-Lx sam- 
ple LF is also much closer we find AM* = 0.64 ± 0.22 and 
Aq = 0.23 ± 0.13 in this case. 

Popesso et al. (2005c) also constructed LFs for galax- 
ies divided into early and late types on the basis of their 
u — r colour. Because of uncertainties in the AT-correction, 
they used a reduced sample of sixteen RASS-SDSS clusters 
located within z < 0.1. As Table 5 shows, Popesso et al. 
(2005c) found a value of M* for the early-type galaxy LF 
intermediate between the values we obtained for our low- 
arid high-Lx cluster samples. In fact, both M* and a for 
this LF are in good agreement with the values we obtain 
for our corresponding high-Lx sample LF. The agreement 
with the low-Lx sample LF is worse, but the M* values 
are still in agreement within < 2a, and the a values differ 
at the ~ 2.3a level. However, for the late-type LF Popesso 
et al. (2005c) find a much brighter characteristic magnitude 
and steeper faint-end slope than we find for both our low- 
and high-Lx cluster samples. We note that unlike the late- 
type LFs produced from 2dFGRS data by Madgwick et al. 
(2002) for the field, and both DP03 and the present work 
for clusters, Popesso et al. (2005c) obtain a brighter char- 
acteristic magnitude for the late-type LF in comparison to 
the early- type LF. The reasons for this discrepancy are not 
clear. 

It is prudent to look for selection effects that could 



2dFGRS: correlation with REFLEX 13 



-21.00 



High-Lj( Sample [19 Clusters] i — ♦- 
Low-L,( Sample [49 Clusters] i o 



-1.40 



: (b) 



8 -1.00 



High-Lx Sample [1 9 Clusters] i — 
Low-Lx Sample [49 Clusters] - ^ 



500 1000 1500 2000 2500 3000 3500 
Acz (km s'^) 



500 1000 1500 2000 2500 3000 3500 
Acz (km s"^) 



Figure 7. The effect on (a) M* and (b) a of limiting selection in redshift space to different slices Acz around each cluster redshift, for 
the luminosity functions of galaxies of any spectral type. The dotted lines show the best-fit value for the appropriate LF parameter and 
sample obtained using the 3 x ctt- selection in cz. 



explain the discrepancy between the early-type luminosity 
functions in our low- and high-Lx cluster samples. In the 
case of the low-Lx sample, which are almost all exclusively 
optically selected, we may worry that they are not in fact 
bound systems, in which case the effect of interloper galaxies 
on the derived LF parameters could be significant. To inves- 
tigate the effect of changing the selection of cluster mem- 
bers in redshift space on the derived LF parameters, we 
constructed LFs using galaxies selected within ±300, ±500, 
±1000, ±2000, and ±3000 km s"^ slices {Acz) around each 
cluster redshift, within the previously calculated selection 
radius of ii2oo (i-e. instead of using 3x the velocity disper- 
sion to select cluster members). Fig. 7 shows the resulting 
Schechtcr function parameters, with la statistical errors de- 
rived from bootstrap resampling, obtained using the same 
method as the results presented in Section 5. We can see 
from this figure that the values of M* and a are not signifi- 
cantly affected by this selection effect, remaining consistent 
with the results presented in Section 5.2. 

We also considered the possibility that a selection ef- 
fect could arise due to galaxy crowding in the centres of the 
high-Lx clusters, because these objects are on average more 
distant than the low-Lx clusters (panel (b) of Fig. 2). There 
are two possible sources of unaccounted-for incompleteness 
in this case: incompleteness in the parent APM galaxy cata- 
logue, due to object blending; and/or incompleteness in the 
2dFGRS spectroscopic catalogue, due to fibre crowding. To 
determine the importance of this effect, we produced lumi- 
nosity functions for the high-Lx sample after excising galax- 
ies within the central 0.25-R200 of each cluster. We present 
the LF for early spectral types in Fig. 8. The derived LF 
parameters (M* = -20.18 ± 0.11, a = -0.86 ± 0.08) are 
clearly consistent with those obtained in Section 5.2, and so 
we conclude that incompleteness effects of this type have no 
significant impact upon our results. 

The physical aperture sampled by a 2dF fibre naturally 
becomes smaller as redshift increases, and so the spectral 
type determined for a particular galaxy may have a mild 
redshift dependence. However, due to seeing effects and fibre 



positioning errors, the effect is likely to be small (see Madg- 
wick et al. 2002). Nevertheless, we attempted to gauge the 
impact of the aperture bias by constructing LFs for both the 
low- and high-Lx samples using only clusters located within 
a common redshift range of 0.03 < z < 0.07. Of the high-Lx 
cluster sample, only 8 clusters of the 19 are found in this 
redshift range, and we found the best fitting early-type LF 
parameters were M* = -19.90 ± 0.13, a = -0.74 ± 0.08. 
The low-Lx sample in this redshift range is made up of 39 
clusters, for which the best fitting Schechter function was 
specified by M* = -20.76 ± 0.13, a = -0.96 ± 0.06. We 
therefore conclude from a comparison of these results with 
those for the full sample that the effect of fibre aperture bias 
is negligible. 

The faint-end of the high-Lx sample is defined by rel- 
atively few clusters. In fact, at magnitudes fainter than 
Mb J > —18.8, only 11 clusters contribute to the LF. Al- 
though the early-type high-Lx sample LF is made up from 
17 per cent more galaxies than the equivalent low-Lx LF, we 
may simply be unfortunate in that the faint-end behaviour 
of these 11 clusters is not typical of the rest of the sam- 
ple. To address this concern, we performed a one-parameter 
fit to the early-type, high-Lx sample LF, with the value of 
the faint-end slope fixed to that found for the low-Lx sam- 
ple, i.e. a = —0.97. The resulting LF is shown in Fig. 9. 
Although the value of M* for the best fit is now brighter 
than found from the two-parameter fit, we find that M* re- 
mains lower than the value obtained from the low-Lx sam- 
ple LF. This result is marginally significant, at the 2a level 
(AM* = 0.30 ± 0.15). Nevertheless, this exercise confirms 
that the low-value of M* found for the high-Lx cluster sam- 
ple arises primarily from the deficit of early-type galaxies 
with Mb J ~ —21 (panel (c) of Fig. 5), rather than from the 
behaviour of the faint-end of the LF. 



6.3 Interpretation 

Since simple selection effects are unable to explain our re- 
sults, we now consider how the difference in the early-type 
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Figure 8. The luminosity function for early (jj ^ —1.4) spectral types for the high-Lx cluster sample, after excising galaxies within 
0.25il2oo of cluster centres. 




Figure 9. Single parameter fit to the luminosity function for early (t] ^ ^1-4) spectral types for the high-Lx cluster sample. The 
faint-end slope has been fixed to a = —0.97, the value found from a fit to the corresponding low-Lx sample LF. 



LF parameters could be driven by environmental differences 

between the low- and high-Lx clusters. Other studies of the 
luminosity function using 2dFGRS data have found evidence 
to suggest that M* brightens and a steepens with increasing 
galaxy density (Croton et al. 2005) or mass (Eke et al. 2004). 
The results presented here seem to suggest that the opposite 
is true when considering the effect of increasing X-ray lumi- 
nosity on the LF. It could therefore be suggested that the 
high-Lx cluster sample is actually composed of low-mass, 
merging systems - Lx is known to be a good tracer of mass, 
but the results of numerical simulations suggest that Lx can 
be boosted during mergers. 

For example, Randall, Sarazin & Kicker (2002, see also 
Ricker & Sarazin 2001) demonstrated that an equal mass 
merger with impact parameter 6 = (i.e. a head-on colli- 
sion) can boost the bolometric X-ray luminosity of the re- 
sultant object by a factor of 5-6. The boost in Lx in this 
case lasts for approximately half a sound crossing time, ~ 1 
Gyr for the high-Lx cluster sample. Without knowledge of 
the merger histories of the individual high-Lx clusters, the 
impact of this effect is difficult to gauge. However, in the 
limiting case of all the high-Lx clusters being boosted in 
Lx by a factor of six, the median unboosted (0.1-2.4 keV) 



X-ray luminosity of the sample is 0.16 x 10*'' erg s~'. By 

comparison, for the majority optically-selected low-Lx clus- 
ter sample the median niaxinmm Lx set by the REFLEX 
flux limit is 0.22 x lO"" erg s"'. Therefore even in this ex- 
treme case both samples would be of similar mass - rather 
than the high-Lx clusters being significantly less massive. 
In addition, as shown in Section 4.2, the median velocity 
dispersion of the high-Lx sample is higher than that for the 
low-Lx sample. These facts imply that the high-Lx clusters 
are more likely to be higher mass systems than the low-Lx 
clusters. 

We suggest that a consistent explanation of the results 
of Sections 4.3 and 5.2 can be achieved by a scenario in 
which the high-Lx clusters are in general older, more dy- 
namically evolved systems than the low-Lx clusters. In this 
case, the higher fraction of early-type, passively evolving 
galaxies present in the high-Lx sample comes about nat- 
urally from there being more crossing times available for 
infalling late-type galaxies to be converted into early-type 
galaxies within the cluster environment. This holds irrc^spec- 
tive of the mechanism responsible for driving the transfor- 
mation, whether it is by interactions with the intracluster 
medium (as in ram-pressure stripping), or by interactions 



between cluster member galaxies (as in the harassment sce- 
nario). 

The excess galaxies that we see at the very bright end of 
the high-Lx cluster LF could be explained by these systems 
having a more extensive merger history than the low-Lx 
clusters, as suggested by the results of Balogh et al. (2002). 
In particular, a 'galactic cannibalism' scenario (e.g. Haus- 
man & Ostriker 1978) - where the brightest cluster galaxies 
grow by the accretion of the successively fainter members 
- could help to bring about the deficit of galaxies seen at 
Mbj ~ —21 in the high-Lx sample, early-type LF, which 
seems to be the cause of the low- value of M* found in this 
case. 

Lin, Mohr & Stanford (2004) studied the it-band LFs of 
93 galaxy clusters divided into two samples of low- and high- 
mass (estimated from the X-ray temperature). Although 
they found the high-mass LF had a brighter characteris- 
tic magnitude than the low-mass LF - in contrast to the 
results of this work for low- and high-Lx clusters - they 
found a deficit of galaxies below M* in the high-mass sam- 
ple in comparison to the low-mass LF. Lin et al. (2004) 
(see also Lin & Mohr 2004) interpreted their result as an 
indication that a major contribution to the masses of the 
very brightest galaxies in clusters comes from mergers with 
other bright (~ M*) cluster galaxies. Perhaps because our 
LFs are constructed purely from spectroscopically-selected 
cluster members - rather than employing a more uncertain 
statistical subtraction of non-cluster members - the effect 
that we see on the derived Schechter function parameters is 
more dramatic. 



7 CONCLUSIONS 

We have correlated the REFLEX catalogue of X-ray selected 
galaxy clusters with the 2dF Galaxy Redshift Survey and 
performed a study of the dependence of cluster galaxy pop- 
ulations upon X-ray luminosity. Using this large, homoge- 
neous combined dataset, supplemented with additional rich 
clusters from the catalogue compiled by De Propris et al. 
(2002), we found: 

(1) An Z/x — Or relation consistent with that derived by 
Ortiz-Gil et al. (2004) for the REFLEX survey team using 
a combination of literature and REFLEX optical follow-up 
data. 

(2) The fraction of early-type, passively evolving galax- 
ies is 12 per cent higher in high-Lx clusters compared to 
low-Lx clusters out to i?2oo, as determined from a volume- 
limited sample of 2dFGRS galaxies with Mbj < — 19, z < 
0.11. In both cases the early-type galaxy fraction falls off 
smoothly with increasing distance from the cluster centre. 

(3) Both the characteristic magnitude M* and faint- 
end slope a of the galaxy luminosity function are lower in 
high-Lx clusters compared to low-Lx clusters. This is caused 
primarily by the low value of M* found for the LF of early 
spectral type galaxies in the high-Lx sample LF, which is 
driven by an underabundance of galaxies with Mj,^ ~ — 21 
relative to the corresponding low-Lx sample LF. 

We believe these results are consistent with a scenario 
in which the high-Lx clusters are more dynamically evolved 
systems than the low-Lx clusters. The higher fraction of pas- 
sively evolving galaxies found in the high-Lx sample could 
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arise through the accretion and conversion of a greater num- 
ber of star-forming field galaxies. The lower value of M* in 
the LF of early spectral types galaxies in high-Lx clusters 
could be explained by 'cannibalism' of bright {Mbj ~ —21) 
members to produce the very brightest cluster galaxies. 
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